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Abstract 



Systems that are not smooth can undergo bifurcations that are forbidden in smooth systems. We review some of the 
phenomena that can occur for piecewise-smooth, continuous maps and flows when a fixed point or an equilibrium 
collides with a surface on which the system is not smooth. Much of our understanding of these cases relies on a 
reduction to piecewise linearity near the border-collision. We also review a number of codimension-two bifurcations 
in which nonlinearity is important. 
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1. Introduction 

A dynamical system is piecewise smooth (PWS) if its phase space can be partitioned into countably many regions 
within which it is smooth. The codimension-one sets on which the dynamics is not smooth are called switching 
manifolds. The recent explosion of interest in these systems is demonstrated by the rapid growth of publications in 
the past two decades shown in Fig. 1, the appearance of a number of texts [1, 2, 3, 4, 5, 6, 7], as well as this issue of 
Physica D. 

Classical bifurcation theory relies heavily on the assumption of a certain degree of smoothness in the singularity 
conditions (e.g. zero eigenvalues), the nondegeneracy requirements (nonvanishing coefficients of some terms in a 
power series), and the use of the center manifold theorem (to achieve a dimension reduction). Bifurcations can occur 
in nonsmooth systems that are forbidden in smooth systems; those that occur due to the interaction of invariant sets 
with a switching manifold are known as discontinuity induced bifurcations. A simple example corresponds to the PWS 
continuous, one-dimensional family of "tent maps". 



for which infinitely many periodic orbits are created as r increases through 1 , without the usual sequence of period 
doubling bifurcations of smooth maps [8, 9].^ 

Of the many types of discontinuities that can occur, perhaps the simplest is a discontinuity in the first derivative; in 
this case the dynamical system is piecewise-smooth and continuous (PWSC). The simplest discontinuity induced bifur- 
cations in these systems correspond to an equilibrium (for ODEs) or a fixed point (for maps) that lies on the switching 
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Figure 1 : Number of papers published each year on piecewise smooth dynamics from a search of the ISI Web of Knowledge® database. Keywords 
were Topic=(piecewise smooth) OR (discontinuity induced) OR (nonsmooth) OR (grazing bifurcation) OR (border-collision bifurcation) AND 
Topic=(dynamics). 



manifold. The resulting bifurcations are known as discontinuous or border-collision bifurcations, respectively. The 
bifurcations occurring in this class of systems are reviewed in this paper. 

More extreme discontinuities include hybrid systems in which the dynamics may be defined by different types 
of models in different regions [11, 12], and Filippov systems where the vector field is discontinuous on a switching 
manifold and the behavior may include sliding along the manifold [13]. 

In this paper we will review some recent results on low-codimension discontinuous and border-collision bifur- 
cations when the switching manifold is smooth at the bifurcation point. One of the difficulties in treating these 
bifurcations generally corresponds to the absence of a center manifold reduction. Even though there are still many 
open questions, there are a number of cases of A^-dimensional flows and maps for which general bifurcation results 
have been obtained. 

We focus on PWSC systems whose components are C* (for some k e N). With this restriction we disregard, for 
instance, maps with square-root terms that arise from regular grazing in Filippov systems [15, 6]; however, we are 
able to exploit a major simplification: border-collision bifurcations can often be treated with a piecewise-linear model. 
For example, nonlinear terms do not influence generic discontinuous saddle-node-like and Hopf-like bifurcations. An 
important consequence is that the size of invariant sets created in these bifurcations grows linearly with parameters 
instead of as some fractional power for the smooth case. 

Piecewise-smooth dynamics arises in many areas including electronics [4], control theory [16], impacting mechan- 
ical systems [5, 17], economics [18, 19, 20] and biology [21, 22, 23] — an extensive set of references can be found in 
[6, 7]. While continuous models are not always appropriate for these applications, there are many cases for which the 
model is given by either a continuous vector field or map. For example, PWSC dynamics occurs in the well-known 
Chua circuit, introduced in 1983, that contains a diode often modeled with a piecewise-linear, current-voltage response 
curve [24]. Similarly power circuitry like DC/DC buck/boost converters, used to decrease or increase voltage levels 
in common devices such as cell phones, can often be modeled by PWSC maps [25]. Piecewise smooth, continuous 
economic models include discrete-time models for trade [26], production-distribution [27] and of a socialist economy 
[28]. Mechanical systems with impacts, or with transitions from static to sliding friction provide another large appli- 
cation area. Continuous models often arise in these systems when a body grazes a rigid surface, leading to so-called 
grazing bifurcations [15, 29]. Biological models that are PWSC include models of neurons [30], neural networks [31], 
and cybernetic models which reduce complex chemical models by selecting optimal reaction pathways [32]. 

In the next section we introduce equations ideal for analyzing local bifurcations of PWSC systems. We solve for 
the equilibria or fixed points for these equations and describe the observer canonical form. Section 3 summarizes recent 
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results regarding codimension-one and two bifurcations for PWSC vector fields. A discussion for the discrete-time 
scenario is presented in §4. 



2. Local Framework 

Near a sufficiently smooth, codimension-one switching manifold one can always choose coordinates so that the 
manifold is represented by the vanishing of one of the coordinates, say the first coordinate [33]. To simplify the 
notation, for any point x e M.^ we let 

s - e[x (2) 

denote the first element of x; consequently, the switching manifold is the plane s = 0. In this case a PWSC dynamical 
system takes different forms in the left-half space, 5 < 0, and the right-half space, s > 0: 

mo-y /(R)(^;^), ,>o • ^^■> 

Here the functions /^'^ : M.^ x M.^ M.'^, i e {L, R) are C'' functions on the A^-dimensional phase space with 
coordinates x and parameters ^ e E*^. We will assume that A: > 1, though often the assumed k will need to be larger 
The function / may represent a vector field, in which case the dynamics is the system of ODEs 

x^fix;^), (4) 

or a map, in which case the dynamics is the discrete-time system 

x>-^ x' ^ f(x;^) . (5) 

An orbit in the discrete case is a sequence {jc, : x,+i - f(x,,^), t e Z). We will discuss these separately in §3 and §4, 
respectively. 

For the simplest bifurcations, we assume that the system has an equilibrium or fixed point on the switching mani- 
fold at some specific parameter value. Without loss of generaUty this point can be taken to be the origin, and we can 
write the parameters as ^ = (fj., rf) where yu e M represents a parameter that vanishes at the border collision. Then we 
can write 

f\x- ij, rf) = ^ib{^l, rf) + Aii^, ?j)x + 0(\xf) + o(k) , (6) 

where Ai and Ar are the N x N Jacobian matrices at the origin. Under the assumption of continuity the two vector 
fields must agree whenever s = 0, which implies that all of the columns of these two matrices must be identical except 
for the first 

Aiei = Ape,- , / 1 . 

The parameter fi corresponds to the primary parameter that unfolds a simple discontinuous or border-collision 
bifurcation, while rj will be used to represent parameters unfolding higher codimension cases. For simplicity, we 
suppress the dependence of the system on j], though the discussion will refer to these additional parameters as needed. 

The nonlinear terms in (6) can be regarded as higher-order in the sense that any structurally stable dynamics of the 
piecewise-linear system will persist when the nonlinear terms are added. When <sc 1, it is useful to consider a scaled 
system, defining z = p when yu and z = x when yu = to obtain 

. ^ f crbiO) + At(Q)z, s<0 
^ \ o-b(O) + An(0)z, s>0 ' 

where cr = sgn(//) (and is when yU = 0). A nice feature of (7) is that fi has been replaced by the discrete parameter 
cr e {-1,0, 1). Note that any structurally stable, bounded invariant set of (7) will approximate an invariant set of (3) 
as yu tends to zero, and moreover that its size (in x) will shrink linearly to zero. This simple feature is one of the most 
prominent characteristics of discontinuous and border-collision bifurcations. 
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If the matrices A,(0) are nonsingular, the imphcit function theorem implies that the ODE (4) has two potential 
equilibria, 

where 

= e] adj(AL(//)) = e] adj(AR(//))) , (9) 

and adj(A) is the adjugate of A (defined by A adj(A) = det(A)/ and extended by continuity for det(A) = 0). Note that 
(8) also pertains to the map (5) if we replace A,- by / - A,- [7]. 

The equilibria (8) are said to be admissible when s*^^^ < and s*^^^ > 0. Otherwise they are virtual. Even virtual 
equilibria play a prominent role in the dynamics of (7), since the eigenspaces of such equilibria typically intersect the 
admissible half-space and these intersections are therefore (forward or backward) invariant sets of the full system. 

A simple consequence of (8) was first noted by Mark Feigin [34] (for the map case): there are two distinct types 
of border collision in this nonsingular case: if det(AL(0)) and det(AR(0)) (or det(/ - Al(0)) and det(/ - Al(0)) for the 
map case) have the same sign, then the equilibria (8) are admissible for different signs of fi, so that near the origin 
the system (3) has a unique equilibrium. On the other hand if the signs of the Jacobian determinants differ, then the 
equilibria are admissible for only one sign of fi and they collide and annihilate as yu — > 0, giving rise to a nonsmooth 
fold bifurcation. 

Additional coordinate transformations can be performed on (7) to reduce the number of parameters. Normally one 
could hope to transform the matrices to Jordan form; however, this cannot always be done simultaneously for both 
matrices, nor would it typically maintain the definition of the switching manifold, s -Q. However, under the condition 
of observability, the matrices can be reduced to companion form 

C = 

where A is a vector and / is the (A^- l)x(A^- 1) identity matrix. Note that the components of A are the coefficients 
of the characteristic polynomial p{A) - det(/l/ - C). Consequently every matrix has a unique companion form. 

Theorem 1 (Observer Canonical Form [6, 7]). There is a coordinate transformation that reduces (3) with (6) to the 
observer canonical form.- (3) with 

f'\x;fi) = iJCN + Ci(n)x + 0{\x\^) + o{k) , (10) 
for small ji and where C,- is the companion matrix ofAj, if and only if Ai{0) has no eigenspace orthogonal to e\. 

Note that the switching manifold for the observer form is still i = 0. The piecewise-linear version of (10) has 
parameters for each matrix, so its complete bifurcation picture is determined by 2N parameters, in addition to 





/ 


-A 








3. Continuous-Time Systems 

The determination of the stability of an equilibrium located on a switching manifold is an important and funda- 
mental yet difficult problem in nonsmooth systems [22, 35]. For PWSC systems, stability depends on the eigenspaces 
of each adjoining smooth subsystem and their relative configuration with respect to switching manifold(s), and, in 
special cases, nonlinear terms of the system [36, 37]. In the absence of nonlinear terms, radial symmetry allows for 
a projection of the dynamics onto a space of dimension one less than the dimension of the system, which is a useful 
simplification [38]. 

When there is a single switching manifold, as for (3), in one or two dimensions the stability problem is straightfor- 
ward [39]. In three or more dimensions the problem is significantly more complicated. For example, an equilibrium 
on a switching manifold may be unstable even if all eigenvalues of both Al and Ar have negative real part [40]. There 
may exist cone-shaped invariant sets akin to slow eigenspaces of smooth systems [41, 42]. By contrast, in two dimen- 
sions such an equilibrium undergoes a Hopf-like bifurcation when, roughly speaking, trajectories of the linearization 



4 



spiral inwards as much as they spiral outwards, per rotation [43, 44, 45]. These considerations are important for some 
models of real-world systems for which some natural constraint results in an equilibrium that remains on a switching 
manifold under parameter perturbation [42, 46]. 

Bifurcations of equilibria in smooth systems of arbitrary dimension are readily classifiable as saddle-node, Hopf, 
and so forth, because the degenerate dynamics occur on a low dimensional subspace (center manifold). For the PWSC 
system, (3), this dimension reduction cannot be performed because the system is not differentiable in a neighborhood 
of the equilibrium. For this reason discontinuous bifurcations resist a simple classification. Though discontinuous 
bifurcations may generate complex dynamics, such as Silnikov homoclinic chaos [47, 48], high-dimensional, PWSC 
systems often can exhibit only low-dimensional dynamics. It remains a difficult problem to develop a general theory 
of dimension reduction for such systems. 

One idea is to look at the eigenvalues of the convex hull of the two Unearizations of the equilibrium [5, 49]: 



For cr e [0, 1], the eigenvalues of A provide a continuous connection between the eigenvalues of Al(0) and Ar(0) and 
are collectively referred to as the eigenvalue path. The idea is that a bifurcation occurs if and only if the eigenvalue 
path intersects the imaginary axis, though currently this remains a conjecture. Furthermore, examples indicate that 
the complexity of the discontinuous bifurcation increases with the number of intersections between the eigenvalue 
path and the imaginary axis. In particular, the presence of exactly one intersection typically yields a discontinuous 
bifurcation that is a natural analogue of a familiar smooth bifurcation. Nonsmooth counterparts of saddle-node, Hopf, 
pitchfork and transcritical bifurcations have been studied [5]. The main difference is that invariant sets created at 
nondegenerate discontinuous bifurcations grow in size linearly, to lowest order. 

When (3) is planar, all possible codimension-one, discontinuous bifurcations that may occur have been determined 
[50, 39]. If the equilibria x**'"* and x*^^\ (8), are nondegenerate they may be classified as saddles, attracting nodes, 
attracting foci, repelling foci or repelling nodes. The bifurcations that result from a border collision between any pair 
of these equilibria are indicated in Fig. 2. 

One possible case is that a collision of an equilibrium with a switching manifold yields no bifurcation in the sense 
that no invariant sets are created at the collision, although there is a piecewise-topological change in the sense of [6]. 
This occurs exactly when the equilibrium fails to change stability at the bifurcation, and is labeled "NB" in Fig. 2. 

When one equilibrium is a saddle and the other is not, the two collide and annihilate in a saddle-node-like bifurca- 
tion, labeled "DSN" in the figure. An additional bifurcation may occur if the second equilibrium is a focus: a periodic 
orbit of stability opposite to the focus can be created at the discontinuous bifurcation [5, 7], this is labeled "DHBSN". 

The discontinuous bifurcation also gives rise to a periodic orbit if one equilibrium is attracting, the other is repelling 
and at least one is a focus, analogously to a Hopf bifurcation and so labeled "DHB". When both are foci the stability 
of the generated orbit is the same as the weaker focus, where strength is measured by the radial factor by which the 
focus takes trajectories toward or away from the equilibrium upon a rotation of 180° [51, 52]. Thus the criticality of 
the bifurcation is determined by linear terms in the ODEs, in contrast to the smooth Hopf bifurcation. 

The case that x*^^^ and x*^^^ are nodes of opposite stability is particularly unusual; we refer to the resulting bifur- 
cation as a "stability-switching" bifurcation, labeled "SS" in the figure. No periodic orbit can be created locally due 
to the presence of the invariant slow manifolds of the nodes. In the absence of nonlinear terms it is useful to add a 
"circle at infinity" to the phase plane, an example is shown in Fig. 3. The resulting discontinuous bifurcation is then 
seen to be analogous to a heteroclinic bifurcation in that at the singularity there exists a continuous loop between the 
equilibrium on the switching manifold and two saddles at infinity. This connection is broken as soon as the equilibrium 
of the piecewise-linear system leaves the switching manifold. The presence of nonlinear terms may cause a periodic 
orbit to be created at the bifurcation; however, this orbit will not shrink to the equilibrium at the singularity. This 
discontinuous bifurcation is seen in, for instance, a piecewise-linear version of the FitzHugh-Nagumo model [53]. 

A variety of codimension-two, discontinuous bifurcations have been analyzed. For example, the solid lines in 
Fig. 2 correspond to codimension-two discontinuous bifurcations for planar systems; their unfoldings are described 
in [7]. Codimension-one discontinuous bifurcations occur along a curve in a two-parameter bifurcation diagram; an 
example taken from [32] is shown in Fig. 4. Curves of classical codimension-one, local bifurcations (such as saddle- 




(1 - o-)Al(0) + o-Ar(O) cr e [0, 1] 
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Figure 2: Codimension-one discontinuous bifurcations involving a single switching manifold for planar systems, taken from [7]. Each row [column] 
represents an equilibrium classification for /''-' [/*'^'], S - saddle; AN - attracting node; AF - attracting focus; RF - repelling focus; RN - repelling 
node. As in [7], names are given to the five distinct codimension-one phenomena: NB - no bifurcation; DSN - discontinuous saddle-node; DHB - 
discontinuous Hopf; DHBSN - discontinuous Hopf-saddle-node; SS - stability switching bifurcation. The chosen layout allows for a classification 
of some codimension-two scenarios — these are the thick lines in the figure. Since smooth parameter variation can produce a change from a saddle 
to a repelling node, the left and right sides of the figure should be considered coincident, as should the top and bottom, see [7] for further details. 



node and Hopf bifurcations) involving an equilibrium of the discontinuous bifurcation may have an endpoint at the 
discontinuous bifurcation curve, like in Fig. 4, since beyond this curve the classical bifurcation is virtual. Roughly 
speaking, as one moves along a curve of saddle-node or Hopf bifurcations, the point in phase space at which the 
bifurcation occurs varies linearly. Equilibria created in a saddle-node bifurcation move apart from one another by an 
amount proportional to the square-root of the parameter change and for this reason the saddle-node bifurcation curve 
is tangent to the discontinuous bifurcation curve at its endpoint [32]. Similarly the amplitude of the periodic orbit 
created in a Hopf bifurcation grows proportional to the square-root of the parameter change, hence this Hopf cycle 
undergoes grazing along a curve tangent to the Hopf bifurcation locus at the codimension-two endpoint [32]. For 
planar systems, at the discontinuous bifurcation on one side of the codimension-two point a periodic orbit is created 
with stability determined by the linear terms of an appropriate piecewise-linear expansion. If this periodic orbit is 
of opposite stability to the Hopf cycle, then the two periodic solutions collide in a saddle-node bifurcation along a 
curve extremely close to the grazing locus; specifically, the series expansions of the two curves agree up to sixth order 
[54]. A rigorous analysis of this bifurcation in higher dimensions, as well as a multitude of other codimension-two 
scenarios, such as the simultaneous occurrence of a discontinuous bifurcation with either a pitchfork or a transcritical 
bifurcation, remains to be done. 
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H<0 n=0 n>0 




Figure 3: Phase portraits including a "circle at infinity" for a planar, piecewise-linear, continuous system undergoing a "stability-switching" 
bifurcation. The equilibria, .v**^' and .r*''^', (8), are nodes of opposite stability and indicated by hollow circles when virtual. When fi > 0, the basin 
of attraction of x'^"^ is bounded by a heteroclinic connection between equilibria at infinity. The trajectory corresponding to the slow eigenvector of 
a:**^* coincides with this connection for x < 0; for jc > this trajectory is virtual and shown by a dotted line. 
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Figure 4: A two-parameter bifurcation diagram of a PWSC model of yeast growth, taken from [32]. Loci of classical saddle-node (SN) and Hopf 
(HB) bifurcations emanate from a locus of discontinuous bifurcations. The Hopf cycle grazes a switching manifold along the dash-dot curve near 
the Hopf locus. 

4. Discrete-Time Systems 

As for continuous-time systems, the stability of fixed points of discrete-time, PWSC systems is typically not 
straightforward to determine. For example the stability of a fixed point at the origin of a planar, piecewise-linear map 
may be treated by reducing it to a circle map with a "dilation ratio" assigned to each point describing the factor by 
which the map takes the point towards or away from the origin [55]. The distribution of the dilation ratios over the 
natural measure of the circle map can reveal the stability of the origin. Of particular interest of late has been the 
scenario referred to as a dangerous bifurcation [56, 57, 58]. This corresponds to the case that the fixed point is stable 
on both sides of the border-collision bifurcation but is not stable at the bifurcation and may be analyzed by considering 
attractors at infinity [59]. 

In one-dimension, the map (3) may be written as 

, ( jj + aix, x<0 
X = <, ^ , (11) 

I yU -H flfiX , X > () 
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and has dynamics that are summarized by Fig. 5, though it should be noted that unstable solutions are not indicated in 
this figure. The map (11) can exhibit discontinuous analogs of saddle-node and Hopf bifurcations for smooth maps, 
the primary diff'erence being that invariant sets separate linearly. Furthermore, (11) can have a stable n-cycle (periodic 
solution of period «) for any «, or a chaotic attractor For further details see [10, 60, 61,6]. 




Figure 5: A two-parameter bifurcation diagram of the one-dimensional, piecewise-linear, continuous map (11), for ^ > 0, taken from [7]. Within 
each shaded region there exists a periodic solution of the indicated period. For < I, there exists a chaotic attractor within the white region. 



Dynamics of piecewise-smooth maps are often easier to describe through the use of symbolics [62, 30, 63, 64]. 
Since here we are considering only dynamics near a single switching manifold, we consider sequences, S, comprised 
of only two letters, L and R. Any orbit, (x,), of (3) with (6), may be assigned a symbol sequence: 



Si 



L , 
R 



if Si < 
if Si > 



where i, = e|jc, and no restriction is placed on Si if = 0. 

For the piecewise-linear case, {xq, xi, . . . x„_i) is an n-cycle if 



where 



Ms ^As„_, ...Aso , 

Ps^ 1 + As„_i + As„_iAs„_, + ■■■+ As„_ 

Following [65, 66, 67], if 1 is not an eigenvalue of the matrix Mg, then 



and 



So 



AstiPs) 



(12) 



(13) 



det(/-M^) 

where - e| adj(/ - Aiijj.)) (as implied by the discussion of (9), which is stated for continuous-time systems). 

These algebraic methods can be used to find admissible periodic solutions, and hence to find regions in parameter 
space where periodic solutions exist for some fixed fi. As in the smooth case, these regions are known as resonance or 
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Arnold tongues, see Fig. 6. Boundaries of resonance tongues correspond to the corresponding «-cycle losing stability 
or undergoing a border-collision. In the latter case one point of the «-cycle typically collides with the switching 
manifold. Remarkably a piecewise-linear expansion of the appropriate n*-iterate of the map will also have the form 
(7), i.e., exactly the same form as the underlying map [60, 7]. Hence (7) may also be used to analyze border-collision 
bifurcations of «-cycles. The stable and unstable manifolds of periodic saddles may undergo homoclinic tangencies 
leading to homoclinic tangles and chaos [68]. 
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Figure 6: Resonance tongues of (7) with (10) for/j > where the multipliers of Cl(0) and Cr(0) are ;-^e*^'"'^L and ^e*-'""", respectively, and 
.vr = 0.8, oji = 0.09, taken from [7]. Period is indicated by the color bar. Within the black regions, forward orbits either approach an orbit of period 
greater than 30 or are aperiodic. 



Resonance tongues in a generic two-parameter bifurcation diagram, like Fig. 6, for piecewise-linear, continuous 
maps often display a structure resembling a string of sausages [69, 70, 71, 18, 68, 65, 72]. A point at which a 
tongue has zero width is known as a shrinking point. When the symbol sequence of the periodic orbit is rotational, 
as defined in [66], the structure near a shrinking point is well understood. Rotational symbol sequences are natural 
as they correspond to rigid translations on an invariant circle. A rigorous unfolding of shrinking point bifurcations 
for rotational periodic orbits can be given under reasonable nondegeneracy conditions. The analysis shows that four 
distinct border-collision curves emanate from a shrinking point; each is tangent to one of the others at the shrinking 
point, an example is shown in Fig. 7. For a nonlinear, PWSC map, shrinking points occur only in the limit fi 0. 
As jj grows, they break apart through the formation of curves of classical saddle-node bifurcations, as shown in Fig. 8 
[67]. The existence of quasiperiodic orbits with irrational rotation numbers can be inferred from the numerics as limits 
of periodic orbits [65]; however, these have not been treated rigorously as of yet, except in some special cases [73, 74]. 

The upper boundary of the period-three tongue (the red region in the lower right corner of Fig. 6) corresponds to a 
loss of stability of a 3-cycle when a multiplier becomes - 1 . At this boundary, a 6-cycle is created, though it is unstable 
and thus not shown in the figure, for further details see [75, 26, 65, 7]. Near the point (wr, ri) ^ (0.366,0.278) small, 
higher-period resonance tongues emanate from a special point on the boundary of the period-three tongue in a manner 
that is only partially understood [65, 7]. 

If one of Al(0) and Ar(0) for (3) has a multiplier on the unit circle, the border-collision bifurcation at fi - Q 
is degenerate. The three codimension-two cases — a multiplier 1, a multiplier -1, or a complex pair of multipliers 
with unit modulus — correspond to the coincidence of a border-collision bifurcation with a saddle-node, a period- 
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0.14 0.16 0.18 0.2 

Figure 7: A magnification of tlie resonance tongue in Fig. 6 corresponding to a rotation number of 1 /7. Tlie boundary of tlie resonance tongue is 
defined by several curves along which there exists a 7-cycle with one point on the switching manifold. 




Figure 8: Resonance tongues for a PWSC map with multipliers and color scheme as in Fig. 6, taken from [67]. Here r/. = 0.2, oj[_ = (jJr, )i = 0.2 
and the map /''"' has a quadratic nonlinearity. The right boundary of each resonance tongue has become C' through the creation of a curve of 
saddle-node bifurcations, and the non-terminating shrinking points have been destroyed. 

doubling, or a Neimark-S acker bifurcation, respectively. Several unfoldings of these bifurcations are given in [7]. 
Recently these three scenarios have been unfolded by Colombo and Dercole [76] without the assumption that the map 
is continuous, in fact without any information about the map on one side of the border-collision bifurcation. The 
results mirror those described in §3 for the continuous-time case. With reasonable nondegeneracy assumptions, loci 
of saddle-node bifurcations and border-collision bifurcations intersect tangentially, whereas loci of period-doubling 
and Neimark-S acker bifurcations intersect the border-collision bifurcation loci transversely. Since the newly created 
solutions along period-doubling and Neimark-S acker bifurcation loci grow in size as the square -root of a parameter, 
they consequently undergo border-collision along a tangent curve as sketched Fig. 9. The nature of the bifurcation at 
this additional border-collision depends on the particular map. For example, it can be shown for the one-dimensional 
case that an attracting chaotic solution may be created along the period-doubling, border-collision curve [77]. 

5. Discussion 

Dynamics related to the interaction of an equilibrium in a PWSC flow, or a fixed point in a PWSC map, with 
a smooth switching manifold are determined by a piecewise series expansion (3) with (6), and, except in special 
cases, only the linear terms of this expansion are important. Consequently invariant sets created at a discontinuous 
or border-colhsion bifurcation grow linearly, to lowest order, with respect to parameter change. Arguably this is 
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Figure 9: A schematic two-parameter bifurcation diagram for a PWSC map. Tlie long daslied curve corresponds to a border-collision bifurcation of 
a fixed point. The solid curves con'espond to loci of classical saddle-node (SN), period-doubling (PD) and Neimark-Sacker (NS) bifurcations. The 
dash-dot curves correspond to the collision of the period-doubled solution and the invariant circle with the switching manifold. The tangency or 
transversality of these curves follows from the general unfolding in [7, 76]. 

the most important feature that distinguishes discontinuity induced bifurcations with bifurcations of smooth systems. 
For example this difference has recently been found in a cardiac model [78, 23] suggesting that a border-collision 
bifurcation and not a classical period-doubling bifurcation generates the observed period-doubled solution. 

Discontinuous bifurcations in one and two-dimensional systems are well understood, as are simple border-collision 
bifurcations. Often the local dynamics is determined by global properties of the appropriate piecewise-linear ap- 
proximation. A number of codimension-two bifurcations that correspond to the coincidence of a discontinuous or 
border-collision bifurcation with one of several classical bifurcations have also been unfolded. However, several 
codimension-two scenarios remain to be fully explored, such as the coincidence of border-collision and Neimark- 
Sacker bifurcations. 

A discontinuity induced bifurcation in a high-dimensional system could be extremely complicated — or it could 
behave like a bifurcation in low-dimensional system. In the latter case, a more complete understanding of the bifurca- 
tion could result from some sort of dimension reduction; however, as of yet no widely applicable theory for dimension 
reduction exists for PWSC systems. 

Resonance tongues in piecewise-linear, continuous maps often have points of zero width (shrinking points) that 
are located at the intersection of several border-collision curves. While much of the structure of shrinking point bifur- 
cations has been explored, there are still a number of open questions. For example, the theory only applies to shrinking 
points in rotational symbol sequences; are there similar bifurcations for non-rotational orbits? What differences arise 
when the lowest order nonlinear terms are not quadratic, but rather some fractional power (as in maps derived from 
sliding bifurcations in Filippov systems)? Also the observation in [65] that curves of resonance tongue shrinking 
points may bound chaotic behavior remains to be explained. Finally, at a shrinking point, the map has an invariant 
polygon that may continue as an invariant circle. This circle can be destroyed by various mechanisms [26, 79, 80]; 
however, analogies of a number of the rigorous results for the smooth case remain to be obtained for PWSC maps, 
perhaps due, in part, to a lack of normal hyperbolicity. 
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